BINARY BLACK HOLE FORMATION AND MERGINGS 



K.A. POSTNOV, M.E. PROKHOROV 

Sternberg Astronomical Institute, 119899 Moscow 
Russia 



The formation and evolution of binary black holes (BH) is studied using the modern evolu- 
tionary scenario for very massive stars. Main sequence stars with masses M > 35Mq are 
assumed to form a BH in the end of their nuclear evolution. The mass of BH formed is 
parametrized as Mbh = kbh x Msn, where Msn is the mass of the pre-supernova star taken 
from evolutionary calculations, kbh < 1- The possibility is explored that a newly formed BH 
acquires a kick velocity 0-250 km/s. Binary BH are found to merge within the Hubble time 
at an appreciable rate only for non-zero kick velocities. We calculate the galactic merging 
rates of binary BH systems, their detection rate by the initial laser interferometers, and the 
distributions of merging binary BH over orbital eccentricities at different frequencies. The 
distribution of angles between BH spins and the orbital angular momentum is also presented. 



1 Introduction 



Among the most promising sources of gravitational waves (GW) which can be detected by the 
initial LIGO/VIRGO/GEO-600 laser interferometers with a sensitivity of h c ~ 10 -21 at the 
frequency 100 Hz, binary black hole (BH) systems play a special role. Their importance stems 
from their larger masses, in comparison to binary neutron star (NS) systems. From observations 
of binary radionulsars and studies of X-ray binaries the average mass of a NS is derived to be 
~ IAMq (see 01), while the mean mass of BH candidates in X-ray binaries is ~ 8M (seei'i). 

Consider a coalescing binary with component masses Mi and Mi lying at a distance r. 
Let the average formation rate per volume of such binaries be T and the merging rate within a 
Hubble time be 1Z, T > 1Z. As the characteristic GW amplitude from such a binary at frequency 
/, which determines the signal-to-noise ratio (<S) at the detector when applying a match filtering 
technique in data analysis, is! h c ~ A1 5/6 ^ _1/6 r _1 (here M = (Mi + M 2 )~ 1/5 (MiM 2 ) 3/5 is the 
"chirp" mass of the binary), the detection rate of such mergings at the detector with a given 
S will be V ~ TZM 15/6 S~ 3 . Thus a strong mass enhancement of binary BH detection rate 
with respect to binary NS detection rate, (Mhhl 'M n s) 15 ' 6 \ cam make binary BH mergings even 



more preferable than binary NS ones U. What is needed is to calculate the formation J- and 
merging rates TL of the corresponding compact binaries. This problem has been addressed earlier 
(e.g. Qa,..., seetl and references therein), but due to the lack of the accurate knowledge about 
BH formation parameters the obtained results have been considered as an order of magnitude 
estimates. 

Binary NS formation and mergings ca n he st udied both using observations of binary ra- 
diopulsars with secondary NS companion or theoretically IISeIb etc. Unlike NS, no 

binary BH system is (and will be) known until GW observatories are at work. 

Difficulties in studying BH formation relate not only to the unknown, from the first princi- 
ples, mass of the star which leaves a BH in the end of its evolution and the mass of the formed 
BH itself, but appears already in the description of the very early evolution of high-mass nor- 
mal stars because of a significant stellar wind mass loss. Until recently, no good evolutionary 
calculations capable of explaining many observational astronomical facts related to massive star 
evolution have been available. The situation, however, seems to have changed after a detailed 
work of Vanbeveren et al. t3. These new calculations of massive star evolution take into ac- 
count the observed stellar wind mass loss from massive stars, and reproduce Wolf-Rayet star 
population characteristics close to those observed in the Galaxy and the Magellanic Clouds. 

Here we present the results of population synthesis of massive binary evolution (see i for 
more detail and references) using these new calculations of Vanbeveren et al. 

2 Galactic binary BH formation, merging, and detection rates 

2.1 A new scenario for massive binary evolution 

A spherically-symmetric mass loss during the main sequence star evolution makes the binary 
semimajor axis a to increase such that after a mass fraction AM has been lost from the system 
a/ao = (Mi + M2)j{M\ + M2 — AM) (e.g. van den HeuveH). The process which can decrease 
the binary separation is the tidal interaction between the components leading to mass transfer 
from one star to another. The most violent mass transfer can result in the common envelope 
stage during which the binary separation can decrease by many times leading sometimes even 
to the coalescence of the components. However, an important fact that follows from the new 
evolutionary calculations El is that primary stars in binaries with initial masses M\ > 4OM can 
not fill their Roche lobes at all because of the strong stellar wind, so it essentially evolves like a 
single star. Of course, in a small fraction of initially very close binaries the more massive primary 
component can fill its Roche lobe and mass transfer can begin. But we assumed that in all initial 
types of binaries the primary with M\ > 40M Q does not fill its Roche lobe. This assumption has 
an advantage of avoiding large uncertainties in the evolution with mass transfer (e.g., the degree 
of non-conservativeness or the common envelope efficiency). As the mass transfer process makes 
the separation between the components decrease, the binary BH merging rates we calculate 
should be considered as lower limits. 

2.2 BH formation parameters 

At present, there is no full understanding of BH formation.Some calculations can be found in 
the literature (see, for example, the recent paper by Fryer til and references therein), but they 
are not self-consistent. Moreover, the high mass-loss pertinent to massive star is ignored in 
these the calculations. So at this stage it seems justified to approximate BH formation by two 
parameters, the initial mass of a main-sequence star leaving a BH in the end of its evolution, 
M cr , and the fraction of the pre-supernova mass that comes to the BH, k^h = M^/Msn- in 
our calculations we have taken M$n from evolutionary calculations 111 and assumed kbh = 0.75 
so that typical BH masses lie within the observed range of BH candidate masses 6-12 Mq. We 



have taken M cr = 35M Q , in accordance with calculations ofEjO. We did not take into account 
a possible BH formation during matter fall-back onto a NS in core-collapse supernpyae (cf. 0) 
or during hypercritical accretion of matter onto NS in a common envelope stage £a. Clearly, 
these additional channels of BH formation will make the formation rate of binary BH systems 
higher by some factor that cannot be relaibly etimated at present. Such accretion-induced BH 
have masses not much higher than the limiting NS mass (~ 3M Q ) and from the point of view 
of GW detection they are closer to NS+BH binaries having smaller chirp masses. 



2. 3 Effect of kick velocity 

One of the most important parameters of binary evolution is the kick velocity acquired by a 
compact star (a NS or BH). The existence of natal kicks during NS formation follows from 
radiopulsar velocity measurements Ii3 and evolutionary studies of pulsar formation (e.g., Ref. 

The physical mechanism for the kick velocity is unknown. The effect of kick velocity on 
the binary compact star formation and coalescence was studied by Lipunov et al. H It is usually 
assumed that the vector of kick velocity is arbitrarily oriented in space and its modulus has 
some distribution f(v) which we used both in a Maxwellian form J{v)m oc v 2 exp[— (v/vq) 2 ] or 
in the form derived from pulsar velocity observations^ J(v)ll ~ (v /t?o)°' 19 (l + (v /t>o) 6 ' 72 ) -1 ^ 2 . 

It is unclear now whether the kick is associated with BH formation as well. However, there 
are astronomical indications that BH candidates (X-ray novae) have higher barycentric velocities 
than ordinary massive binaries 110, which can be explained either assuming a substantial mass 
loss (which we do not think to occur) or by a kick velocity of ~ 100 km/s during BH formation 
El. A long-term periodicity observed in some BH candidates can be also due to spins of the 
component being non-parallel to the orbital angular momentum. So in our calculations we 
assumed BH kick velocities distributed in the same way as for NS but with Vq 1 = Vq s (1 — 
kbh)/(l — Mov/Msn), where kbh = Mbh/MsN is the fraction of the pre-supernova mass M$n 
collapsing into BH, Mov = 2.5Mq is the upper mass limit of NS. Clearly, for kbh = 1 no 
additional velocity is acquired by BH, and when Mbh —> Mov the kick velocity Vbh v ns . 

The effect of kick velocity is twofold. Depending on the value and the kick orientation 
relative to the velocity of the pre-supernova star before the explosion, it can either disrupt or 
bind stronger the system (see e.g. Yamaoka et al.il for more detail). We found that for binary 
BH with typical masses of 7 — 10M Q which form during the evolution of massive binaries the 
kick velocity is mostly affects the merging rate. Without kick velocity double BH form butmost 
of them recide in too wide binaries to merge in a time shorter than the Hubble time (cf. ca), so 
binary BH merging rate 1Z is very small to be of significance for GW observations. 

The situation, however, strongly changes for non-zero kicks during BH formation. Then a 
part of the double BH systems are formed vith very high eccentricities e ~ 1 after the second 
BH formation. As is well known Ea, the eccentricty shortens the coalescence time due to GW 
emission t c (eo,M) of a binary with the chirp mass M, initial orbital period Pq, and orbital 
eccentricity eo 

t c (e ,M) =t (M,Po)He Q ), (1) 
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rapidly tends to zero when e — ► 1. A high eccentricity thus enables binary BH systems with 
very large orbital semiaxes (periods) to merge in a Hubble time. If the kick velocity is higher 
than the orbital velocity of the exploding component, it preferentially makes the binary system 
unbind, so the dependence of binary BH merging rates on the kick velocity is expected to have 
a maximum. This effect is confirmed by Monte-Carlo modeling (see Fig. [l]). 

3 Population synthesis modeling of BH+BH mergings 

The results of population synthesis calculations of binary BH merging rate in a model galaxy of 
1O 11 M0 (assuming all stars forming in binaries) with a constant star fromation rate are shown 
in Fig. 1. The binary evolution parameters were taken in the conventional form: Salpeter initial 
mass function for the primary mass Mi, a flat mass ratio q = M2/M1 distribution, a flat initial 
semiaxes distribution a (seeH for more detail). The galactic binary BH merging rates are shown 
as a function of the kick velocity during BH formation. The calculations were done with a delta- 
function like kick velocity distribution (the use of more complex distributions do not change the 
results significantly and not shown in the Figure). In Fig. 1 (right panel) the detection rate 
of binary BH coalescences by the initial laser interferometers (h c = 10 -21 at u = 100 Hz) in 
one- year integration is shown as a function of BH kick velocity. It is seen that both the galactic 
merging rate and detection rate of binary BH systems rapidly increases with the assumed kick 
velocity amplitude and reach a maximum of 1Z ~ 2.5 x 10~ 5 yr _1 and T> ~ 20 detections per 
year, correspondingly, at vt — 120 km/s. Since V ~ M 15/6 K, the TZ(v k ) and V(v k ) dependences 
have similar shapes. 

In Fig. 2 we present the distribution of eccentricities of merging BH+BH systems at the 
characteristic frequencies few = 10 -2 Hz and vqw = 10 Hz for LISA and LIGO/VIRGO/GEO- 
600 interferometers. The eccentricity distributions has a power-law shape /(e) ~ e~ 13 in a wide 
range, which reflects the transformation of the initial eccentricity in the course of the orbital 
evolution of a binary system of two point masses due to gravitational wave emission and the 
distribution of the intial eccentricities (i.e. immediately after the second BH formation) of 
the merging binaries. For example, for a flat initial eccentricity distribution we would obtain 
/(e) ~ e - ^ 2 / 19 , which directly follows from the relation between semimajor axis and eccentricity 
evolution o. A small fraction of merging binary BH do has appreciable eccentricities. These 
systems evolve from wide, extremely eccentric binaries. 

Unlike a spherically symmetric explosion of one of the components in a binary system, 
when a sudden mass loss AM influences only on the binary major semiaxis and eccentricity, the 
asymmetric explosion with kick velocity also changes the space orientation of the orbital angular 
momentum. Unless the vector of the kick is non-central, the spin axis of the star is not changed. 
Initially, each component of a binary should have spins parallel to the orbital angular momentum. 
The spins are assumed to remain parallel to each other during the evolution. After each BH 
formation with kick the orbital angular momentum changes its direction, and the resulting 
distribution of cosines of the angle between the BH spins and the orbital angular momentum 
(denoted by cos J) after the second BH formation is presented in Fig. 3. Remarkably, even for 
small kicks of a few ten km/s an appreciable fraction (30-50%) of the merging binary BH should 
have cos J < 0. Note also a non-monotonic change of the distribution form with kick. Only a 
tiny fraction of binaries can have spins antiparallel to the orbital angular momentum. 

4 Conclusions 

We have carried out calculations of binary BH merging rate and detection rate by the initial 
laser .GW interferometers in the framework of the new evolutionary scenario for very massive 
stars E-3 . We have found that both the merging and detection BH+BH rates are significant 




Figure 1: Left: BH+BH merging rates calculated for a 1O 11 M galaxy with a constant star formation rate, as a 
function of assumed kick velocity during BH formation with M cr = 35M Q , for kbh = 0.5 and 0.75. Right: Detection 
rate per year of BH+BH mergings by the initial LIGO/VIRGO/GEO-600 laser interferometers (h c = 10~ 21 at 
/ = 100 Hz), as a function of kick velocity during BH formation. 




eccentricity eccentricity 

Figure 2: Differential distribution of the orbital eccentricities of merging BH+BH systems at frequency 10~ 2 Hz 
(left panel) and 10 Hz (right panel). Figures mark the kick velocity amplitude. 




Figure 3: Integral distribution of the angles between the orbital angular momentum and BH spins in merging 
BH+BH systems. Figures mark the kick velocity amplitude. 



from the point of view of GW detection by the initial LIGO/VIRGO/GEO-600 detectors only if 
non-zero kick velocity is assumed during BH formation. This important point and the degree of 
compatibility of the new scenario with existing astronomical observations surely deserves further 
studies, which are in progressEl. We confirm the previous result iH that binary BH still remain 
the most promising sources for the initial laser GW interferometers. 

Our calculations show that in this scenario, merging binary BH can have non-zero eccentric- 
ities at frequencies 1(T 2 Hz and 10 Hz relevant to LISA and LIGO/VIRGO/GEO-600 frequency 
bands. We also found that even small kick velocities during BH formation result in an apprecia- 
ble angle between the BH spin axes and the orbital angular momentum for a significant (30-50%) 
fraction of coalescing binary BH. These findings should be taken into account in constructing 
GW templates for data analysis from GW detectors. 
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